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g ; Abstract 

O 

We explore a version of the cosmological dilaton-fixing and decoupling mech- 
anism in which the dilaton-dependence of the low-energy effective action is 

Q-l' 9 rh 

extremized for infinitely large values of the bare string coupling g s = e' p . We 



X 



study the efficiency with which the dilaton cf) runs away towards its "fixed 
point" at infinity during a primordial inflationary stage, and thereby approx- 
imately decouples from matter. The residual dilaton couplings are found to 
be related to the amplitude of the density fluctuations generated during infla- 
tion. For the simplest inflationary potential, V(x) = ^rn^(4>)x 2 : the residual 
dilaton couplings are shown to predict violations of the universality of gravita- 
tional acceleration near the Aa/a ~ 10 -12 level. This suggests that a modest 
improvement in the precision of equivalence principle tests might be able to 
detect the effect of such a runaway dilaton. Under some assumptions about 
the coupling of the dilaton to dark matter and/or dark energy, the expected 
time-variation of natural "constants" (in particular of the fine-structure con- 
stant) might also be large enough to be within reach of improved experimental 
or observational data. 
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I. INTRODUCTION 



All string theory models predict the existence of a scalar partner of the spin 2 gravi- 
ton: the dilaton 0, whose vacuum expectation value (VEV) determines the string coupling 
constant g s = e 9 ^ 2 [1]. At tree level, the dilaton is massless and has gravitational-strength 
couplings to matter which violate the equivalence principle [2]. This is in violent conflict 
with present experimental tests of general relativity. It is generally assumed that this con- 
flict is avoided because, after supersymmetry breaking, the dilaton might acquire a (large 
enough) mass (say > 10~ 3 eV so that observable deviations from Einstein's gravity are 
quenched on distances larger than a fraction of a millimeter). However, Ref. [3] (see also [4]) 
has proposed a mechanism which can naturally reconcile a massless dilaton with existing 
experimental data. The basic idea of Ref. [3] was to exploit the string-loop modifications of 
the (four dimensional) effective low-energy action (we use the signature — h ++) 

S = J Jxy/jj R+^M- [2D0 - (V0) 2 ] - \ B F {<f>) F 2 - V + • • •) , (1.1) 

i.e. the 0-dependence of the various coefficients -Bj(0), i = g, 0, F, . . . , given in the weak- 
coupling region (e^ — > 0) by series of the form 

5,(0) = e~* + cf + cf + cf e 2 ^ + • • • , (1.2) 

coming from genus expansion of string theory: Bi = S n g^^c^, with n — 0, 1, 2, It was 

shown in [3] that, if there exists a special value m of <fi which extremizes all the (relevant) 
coupling functions i^ rl (0), the cosmological evolution of the graviton-dilaton- matter system 
naturally drives towards m . This provides a mechanism for fixing a massless dilaton at 
a value where it decouples from matter ("Least Coupling Principle"). A simple situation 
where the existence of a universally extremizing dilaton value <p m is guaranteed is that of S 
duality, i.e. a symmetry g s <-> l/g s , or — > —0 (so that m = 0). 

It has been recently suggested [5] that the infinite-bare-coupling limit g s — > oo (0 — > +oo) 
might yield smooth finite limits for all the coupling functions, namely 

B i ((f>)=C i + 0(e-1 > ). (1.3) 

Under this assumption, the coupling functions are all extremized at infinity, i.e. m = +oo. 
The late-time cosmology of models satisfying (1.3) has been recently explored [6]. In the 
"large N"-type toy model of [5] it would be natural to expect that the 0(e^^) term in 
Eq. (1.3) be positive, so that Bi((f>) be minimized at infinity. This would correspond to 
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couplings Xi((f>) ~ = C^ 1 — 0(e~^) which are maximized at infinity. Note, however, 

that the most relevant cosmological coupling for this work, the coupling to the inflaton, A(0), 
contained in V (see Eq. (2.11) below) is closer to a Bi than to its inverse. Thus \(<t>) is 
naturally minimized at infinity (see further discussion of this point below), a crucial property 
for the attractor mechanism of [3,4]. 

In this paper* we shall consider in detail the early-time cosmology of models satisfy- 
ing (1.3). More precisely, our main aims will be: (i) to study the efficiency with which 
a primordial inflationary stage drives towards the "fixed point" at infinity m = +oo 
(thereby generalizing the work [8] which considered the inflationary attraction towards a 
local extremum <f> m ), and (ii) to give quantitative estimates of the present violations of the 
equivalence principle (non universality of free fall, and variation of "constants"). Our most 
important conclusion is that the runaway of the dilaton towards strong-coupling (under the 
assumption (1.3)) naturally leads to equivalence-principle violations which are rather large, 
in the sense of not being much smaller than the presently tested level ~ 10~ 12 . This gives 
additional motivation for the currently planned improved tests of the universality of free fall. 
Within our scenario, most of the other deviations from general relativity ( "post-Einsteinian" 
effects in gravitationally interacting systems: solar system, binary pulsars, ...) are too small 
to be of phenomeno logical interest. However, under some assumptions about the coupling 
of to dark matter and/or dark energy, the time variation of the natural "constants" (no- 
tably the fine-structure constant) predicted by our scenario might be large enough to be 
within reach of improved experimental and/or observational data. The phenomenologically 
interesting conclusion that equivalence-principle violations are generically predicted to be 
rather large after inflation (in sharp contrast with the results of [8]) is due to the fact that 
the attraction towards an extremum at infinity is much less effective than the attraction to- 
wards a (finite) local extremum as originally contemplated in [3] . This reduced effectiveness 
was already pointed out in Ref. [4] within the context of equivalence-principle-respecting 
tensor-scalar theories (a la Jordan-Fierz-Brans-Dicke). 

II. DILATON RUNAWAY 

In this section we study the dilaton's runaway during the various stages of cosmological 
evolution. We first show (subsection II A) that, like in the case of a local extremum [8], 
inflation is particularly efficient in pushing towards the fixed point. We will then argue 



*The main results of this work have been recently summarized in a short note [7]. 
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(subsection II B) that the order of magnitude of the bare string coupling ~ e C(p does not 
suffer further appreciable changes during all the subsequent evolution. 



A. The inflationary period 



Assuming some primordial inflationary stage driven by the potential energy of an inflaton 
field x, and taking into account generic couplings to the dilaton 0, we consider an effective 
action of the form 



= / d^g [250 - (V0) 2 ] - \ 5 x (0)(Vx) 2 - V(x, <t>) ) ■ (2.1) 



In this string-frame action, the dilaton dependence of all the functions -Bi(0), V(x, 4>) is 
assumed to be of the form (1.2). It is convenient to replace the (a-model) string metric g^ v 
by the conformally related Einstein metric g^ u = C B g {(f))g^ ui and the dilaton field by the 
variable 



dd> 



4 \bJ B g 2 B g 



, B' = dB/d<f> . 



(2.2) 



The normalization constant C is chosen so that the string units coincide with the Einstein 
units when — > +oo: C B g (+oo) = 1. [Note that C — 1/C g in terms of the general notation 
of Eq. (1.3).] Introducing the (modified) Planck mass 



m 2 P 



1 



4nG Ca> ' 



(2.3) 



and replacing also the inflaton by the dimensionless variable \ — C rrip 1 x, we end up 
with an action of the form 



= J d A x^g 



^fR-^f (Vv,) 2 - ^ FM(Vx) 2 - m% V( X , <p) 



(2.4) 



where 



F(<p) = B x {4>)/B g {4>) , VixM = C- 2 mp 4 R- 2 (0) V(x,<f>) 



(2.5) 



In view of our basic assumption (1.3), note that, in the strong-coupling limit — > +oo, dtp/ 'd<ft 
tends, according to Eq. (2.2), to the constant (C < j } /2C g ) 1 ^ 2 , while the dilaton-dependent factor 
F(tp) in front of the inflaton kinetic term tends to the constant C x /C g . The toy model of 
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Ref. [5] suggests that the various (positive) constants Cj in Eq. (1.3) are all largish and 
comparable to each other. We shall therefore assume that the various ratios Ci/Cj are of 
order unity. The most important such ratio for the following is c = (2C g /C ( p) 1 ^ 2 which gives 
the asymptotic behaviour of the bare string coupling as 

9 2 s =e*~ e c * . (2.6) 

In view of the fact that, in the strong-coupling limit we are interested in, the factor F(ip) in 
Eq. (2.4) quickly tends to a constant, we can simplify our analysis (without modifying the 
essential physics) by replacing it by a constant (which can then be absorbed in a redefinition 
of x)- Henceforth, we shall simply take F{tp) = 1. [See, however, the comments below 
concerning the self-regenerating inflationary regime.] 

Following [4,9] it is then useful to combine the Friedmann equations for the scale factor 
a(t) during inflation (ds 2 = —dt 2 + a 2 (t) 5ij dx l dx j ) with the equations of motion of the two 
scalar fields x(*)> V?(0> to write an autonomous equation describing the evolution of the two 
scalars in terms of the parameter 

p = J Hdt = J ^ dt = In a + const (2.7) 

measuring the number of e-folds of the expansion. For any multiplet of scalar fields, p = (y? a ), 
this yields the simple equation [4,9] 

_J_ v " + 2 p' = -Vp In \V(<p) | , (2.8) 



3 — ip 

where ip' = dcp/dp, and where all operations on <p are covariantly defined in terms of the 
cr-model metric defining the scalar kinetic terms [da 2 = 'Jabif) df a df b ). In our simple model 
(with F((p) = 1), we have a flat metric da 2 = dip 2 + dx 2 - [Note that, when 7a&(y?) is curved 
the acceleration term <p" involves a covariant derivative.] 

The generic solution of Eq. (2.8) is easily grasped if one interprets it as a mechanical 
model: a particle with position (p, and velocity-dependent mass m(ip') = 2/(3 — y' 2 ), moves, 
in the "time" p = lna + est, in the manifold da 2 under the influence of an external potential 
In an d a constant friction force —lip'. If the curvature of the effective potential 

In |V(^)| is sufficiently small the motion of (p rapidly becomes slow and friction-dominated: 

2^ ~ -VinlnVtcp). (2.9) 
dp 

Eq. (2.9) is equivalent to the usual "slow roll" approximation. 
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Consistently with our general assumption (1.3), we consider potentials allowing a strong- 
coupling expansion of the form: 

VixM = V (x) + KixK^ + 0(e~ 2c n , (2.10) 

where V (x) is a typical chaotic-inflation potential with V^(0) = 0, while 14(0) = v i > can 
possibly provide (if v\ > 0) the effective cosmological constant driving today's acceleration 
in the scenario of [6]. For the sake of simplicity we shall discuss mainly the "factorized" 
power-law case V (x) ~ Vi(x) ~ X n f° r which we can conveniently write V in the form: 

n 

vixM = K<p) — > (2-n) 

n 

with a dilaton-dependent coupling constant X((p) of the form 

A( V 9) = A 0O (l + 6 A e^). (2.12) 

This example belongs to the class of the two-field inflationary potentials discussed in [10]. We 
have checked that our results remain qualitatively the same for the more general potential 
(2.10) provided that Vo(x) and V\{x) are not extremely diffent and given the fact that v\ 
is phenomeno logically constrained to be very small. [Note that, within the simplified model 
(2.12), the ratio V 1 {x)/Vq{x) is equal to the constant coefficient b\] 

The universal (positive) constant c appearing in the exponential e~ cip is the same as in 
Eq. (2.6) [i.e. c = (2C g /C^ ) ) 1 ^ 2 , which is expected to be of order unity]. The coefficient 
b\ in (2.12) is such that b x eT Cip ^ b\e~^ roughly corresponds to a combination of terms 
~ ± Ci l 0(e~^) coming from the strong-coupling asymptotics of several -B«(0), Eq. (1.3) (see 
Eq. (2.5)). In the toy model of [5] one would therefore expect 6a to be smallish. Anyway, we 
shall see that in final results only the ratios of such fej coefficients enter. More important than 
the magnitude of b\ is its sign. It is crucial for the present strong-coupling attractor scenario 
to assume that b\ > 0, i.e that \(<p) reaches a minimum at strong-coupling, <p — > +oo. 
Note again that this behaviour is consistent with the simple "large iV"-type idea of [5] if we 
assimilate X((p) to one of the inverse couplings Bi appearing in (1.1) (for instance Bp ~ gp 2 , 
where gp is a gauge coupling), rather than to the coupling itself. If the latter were the case, 
\(tp) would reach a maximum as <fi — > +oo, and the attractor mechanism of [3] would drive 
towards weak coupling (<f) — > — oo). However, the Einstein-frame ^-dependence of V(x) 
gets contributions from several Bf n ((f)), Eq. (2.5), which might conspire to minimize it at 
strong coupling. This feature is also probably necessary in order to solve the cosmological- 
constant problem through some argument by which the vacuum at infinity has vanishing 
energy density. 
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Substituting the potential (2.10) into the slow roll equation (2.9) and assuming (for 
simplicity) that Vi(x)e~ c<p is significantly smaller than V (x) leads to a decoupled set of 
evolution equations for x an d V 9 (where V = dV/dx)' 

**-- X -YL (213) 
dp~ 2 14' [ } 

^ = i ce -^. (2.14) 
rfp 2 Vb 

Given some "initial" conditions Xin, <Pm (discussed below) at some starting point, say p — 0, 
the solution of Eqs. (2.13), (2.14) is simply 



e^-e c ^+ - /^ Vl(x(p)) (2 16) 



which simply become: 

P = 1 (X?n - X 2 ) , e^ + b -^ X 2 = const. = + ^ x? n , (2.17) 



in the simplified case of eqs. (2.11), (2.12). 

Equations (2.17) show that, in order for the string coupling g 2 s ~ e CLp to have reached 
large values at the end of inflation, a large total number of e-folds must have occurred while 
the (dimensionless) inflaton field x decreases from a large initial value, to a value of order 
unity (in Planck units). To get a quantitative estimate of the string coupling at the end of 
inflation we need to choose the initial conditions Xin, 'Pin- A physically reasonable way (which 
is further discussed below) of choosing Xin is to start the classical evolution (2.13)-(2.17) at 
the exit of the era of self- regenerating inflation (see [10] and references therein). We will now 
show how to relate the exit from self-regenerating inflation to the size of density fluctuations 
generated by inflation. 

Let us recall (see [10] and references therein) that the density fluctuation 5 = Sp/p 
on large scales (estimated in the one-field approximation where the inflaton x is the main 
contributor) is obtained by evaluating the expression 



2 



4 1 /2\5 V 3/ ' 
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at the value x = Xx, at which the physical scale we are considering crossed the horizon 
outwards during inflation. For the scale corresponding to our present horizon this usually 
corresponds to a value Xx(Ho) {xh for short) reached some 60 e-folds before the end of slow- 
roll. From Ref. [10], xh — 5y/n for the model (2.11) (and with our modified definition of x)- 
The numerical value of Sh = S(xh) which is compatible with cosmological data (structure 
formation and cosmic microwave background) is Sh — 5 x 10~ 5 . In the model (2.11) the 
function S(x) defined by (2.18) scales with x a s X~ ■ Putting together this information 
we obtain a relation between Xm an d S(xm), which involves the value of the observable 
horizon-size fluctuations Sh = S(xh)' 



<KXin) _ (x 

Kxh) \Xh, 



n+2 
2 

in \ 



(2.19) 



i.e. 





where we introduced the short-hand notation S m = S(xm)- 

Inserting Eq. (2.20) into Eq. (2.16) we then obtain the following estimate of the string 
coupling constant after inflation as a function of y? in and S(xm) 

r 2 r 2 25r 2 ( S- \ 

e c cnd _ e c in „ - {Vl /V )p ~ ^(Vi/V )xl - — (K/Vo) i^j , (2.21) 

where (Vi/Vo) denotes the average value of V\/Vq\ (Vi/Vq) = / dp(Vi/Vo)/ J dp [note that 
this average ratio is equal to b\ in the simplified model (2.12)]. 

To get a quantitative estimate of e c<Amd we still need to estimate the value of S(xin) 
corresponding to the chosen "initial" value of the inflaton. As we will now check, taking for 
Xm the value corresponding to the exit from self-regenerating inflation corresponds simply 
to taking S(x i n ) ~ 1. Indeed, let us first recall that, during inflation, each (canonically 
normalized) scalar field (of mass smaller than the expansion rate H) undergoes typical 
quantum fluctuations of order H/(2tt), per Hubble time (see, e.g., [10]). This implies (for 
our dimensionless fields) that the value of x &t the exit from self- regeneration, say Xex, is 
characterized by H ex / (2ir) [d x V/(2V)] cx , where H = H/rhp is the dimensionless Hubble 
expansion rate and where the right-hand side (RHS) is the classical change of x P er Hubble 
time (corresponding to the RHS of Eq.(2.13)). Using Friedmann's equation (in the slow- roll 
approximation) H 2 X (2/3)V(x C x), h is easily seen that that the exit from self-regeneration 
corresponds to <5(Xcx) ~ 4/3 ~ 1. It is, a posteriori, physically quite reasonable to start using 
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the classical evolution system only when the (formal extrapolation) of the density fluctuation 
5(x) becomes smaller than one. 

Within some approximation (see [10]), one can implement the effect of the combined 
quantum fluctuations of (</?, %) by adding random terms with r.m.s. values H/2n on the 
right hand side of equations (2.13) and (2.14), dx/dp and dip /dp being precisely the shifts 
of the fields in a Hubble time. The system of equations becomes thus of the Langevin-type 

where H [(2/3) V(x, 0)] 1//2 (in the slow- roll approximation) is the dimensionless expansion 
rate, and where £i and £2 are (independent) normalized random white noises: 

(UPi)UP2)) = S ij 8(p 1 -p 2 ), i,j = 1,2. (2.24) 



When the random force terms dominate the evolution in either equation (2.22) or equation 
(2.23) the quasi-classical description (2.13) , (2.14) breaks down. The phase space of the 
sytem can thus be roughly divided into four regions according to whether the evolution of 
none, one or both of the two fields is dominated by quantum fluctuations. This is depicted in 
Figure 1 where such regions are delimited by dashed, thick curves in the case of a power-law 
potential (2.11). 

Apart from factors of order one, the evolution of the inflaton x is quasi-classical in the 
region under the line x — ^^ n+2 ' > - 111 the chaotic inflationary models [10] such an inflaton's 
value corresponds to the exit from the self-regenerating regime and to the beginning of the 
quasi-classical slow-roll inflation. As mentioned above it also corresponds to a perturbation 
S(x) ~ 1- Somewhat surprisingly, in the model at hand, however, the quasi-classical region 
for the inflaton evolution x ~ ^«) 1 ^™ +2 ' ) is affected by quantum fluctuations that still dom- 
inate the evolution of ip in the region above the hyperbola-like curve x = b 2 / 71 \^/ n e~ 2c<f / n . 
We must therefore study, in some detail, the evolution of the system in the presence of the 
noise term for (p as in eq. (2.23), assuming a quasi classical evolution for x- This is done 
in the appendix. The final result is that, if we start at x ~ ^^^ n+2 \ (classical evolution 
for x) the average value of e C(p is multiplicatively renormalized, by a factor of order unity, 
with respect to the classical trajectory e c<fcl , given by solving equations (2.13) - (2.14), i.e. 
(e C(p ) = 0(l)e ciPcl . One also finds that the dispersion of e cip around its average value is 
comparable to its average value. 
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50 100 150 200 250 300 

exp(c(/>) 

FIG. 1. The phase space of the system is represented in the case of a power-law potential 
(2.11) with n = 2, b x = 0.1 and = 1(T 10 . The thick-dashed (red) curves delimitate the 
quantum behaviour of the two helds, the horizontal curve x = Aoo and the hyperbola-like 

curve x = b^ n \oo^ n e~ 2cip / n being the limit of the quantum behaviour for x and (p respectively. In 
the white region both helds have a classical behaviour. The last "fully classical" trajectory has been 
represented by a thick (blue) curve. The bright-gray regions are those where either the <p or the 
X evolution are dominated by quantum fluctuations. The fully-quantum region is the dark-gray 
region on the top right. 

We shall not try to discuss here what happens in the self-regenerating region \ > Xm ~ 
A ( ^ 1 /( n+2 ' ) . Let us recall that the simple decoupled system (2.13) , (2.14) was obtained by 
neglecting the kinetic coupling term F(ip) in Eq. (2.4). If we were to consider a more general 
model, we would have more coupling between x an d y? and we would expect that (contrary 
to Fig. 1 which exhibits a "classical ip region" above the "quantum x line" ) the evolution in 
the self-regenerating region involves a strongly coupled system of Langevin equations. Then, 
as discussed in [10], solving such a system necessitates to give boundary conditions on all 
the boundaries of the problem: notably for x ~ > °°> but also for ip — > +oo and ip — > — oo. 
We leave to future work such an investigation (and a discussion of what are reasonable 
boundary conditions). In this work we shall content ourselves with "starting" the evolution 
on the quantum x boundary line Xm with some value ip = (pi n , assuming that e cipin is smaller 
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than the driving effect due to inflation, i.e. than the RHS of Eq. (2.21). [This assumption is 
most natural in a work aimed at studying the "attracting" effect due to primordial inflation.] 

Going back to our result (2.21), we can now insert, according to the preceding discussion, 
the values 5- m = 1 and e cipin <C e CiPcnd . Finally, in the simplified model (2.11), (2.12), we get 
the estimate: 

e^cnd = . e ^ c ,,end _ . bx (§ H y—2 . (2.25) 

A more general analysis, based on the potential (2.10) leads to the same final result but 
with n replaced by some average value of x ^°'* , and with b\ replaced by some average of the 
ratio Vi/Vq. Note that smaller values of the exponent n lead to larger values of e ciPcnd , i.e. to 
a more effective attraction towards the "fixed point at infinity" . The same is true if we take 
different exponents Uq and n\ (for Vq and V\ respectively) and assume Vi(xi n ) ^> Vo(Xin) to 
hold as a result of Xm ^> 1 and n\ > uq. Also note that, numerically, if we consider n = 2, 
i.e. the simplest chaotic-inflation potential V = | ra^) x 2 > Eq. (2.25) involves the large 
number 12.5 x 5jj l ~ 2.5 x 10 5 . In the case where n = 4, i.e. V = \ A(v?) x 4 , we have instead 
the number 12.5 x d^ 2 ^ 3 ~ 0.92 x 10 4 . To understand the phenomenological meaning of these 
numbers we need to relate e c<Amd to the present, observable deviations from general relativity. 
This issue is addressed in Section III after having argued that the post- inflationary evolution 
of ip is sub- dominant. 



B. Attraction of ip by the subsequent cosmological evolution 

We have discussed above the efficiency with which inflation drives the dilaton towards 
a fixed point at infinity. We need to complete this discussion by estimating the effect of 
the many e-folds of expansion that took place between the end of inflation and the present 
time. To address this question, we need to study in more detail the coupling of a runaway 
dilaton to various types of matter, say a multi-component distribution of (relativistic or 
non-relativistic) particles. We work in the Einstein frame, with an action of the type 



(2.26) 

~J2 m A [(p(x A )]^-g fJ , l/ (x A )dx^dx A . 

A J 

Following [2,3], one introduces the crucial dimensionless quantity 
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d In rriA(<p) 
dip 



(2.27) 



measuring the coupling of <p to a particle of type A. [For consistency with previous work, we 
keep the notation a a but warn the reader that this should not be confused with the various 
gauge coupling constants, often denoted ctj = g 2 /4n.} The quantity a a determines the effect 
of cosmological matter on the evolution of p through the general equation [4,3] 

^ + ( 1 - - V = - £ M<f) , (2-28) 

where the primes denote derivatives with respect to p = In a + const and where p = Ha pa 
and P = Ha Pa are the total "material" energy density and pressure respectively, both 
obtained as sums over the various components filling the universe at the exception of the 
kinetic energy density and pressure of tp, p k = (m 2 P /2)(dp/dt) 2 = im 2 p /2)H 2 p' 2 and P k = p k . 
Accordingly, the Friedmann equation reads 

3H 2 = 4-Ptot = ^r + H 2 p' 2 . (2.29) 
mp mp 

Note that p and P may also account for the potential energy density and pressure of 
the scalar field, py = V(p), Py = —py and that one can formally extend Eq. (2.28) to 
the "vacuum energy" component V(ip) by associating to the potential V(ip) the mass scale 
m v (<p) = V{ip)^ which gives ay = d\nmy(ip)/dip = \ d\nV(p)/dip. Equation (2.8) is then 
recovered in the limit where the scalar field is the dominant component. 

In the simple cases (which are quite frequent; at least as approximate cases) where one 
"matter" component, with known "equation of state" Pa/ Pa = wa = const., dominates 
the cosmological density and pressure, Eq. (2.28) yields an autonomous equation for the 
evolution (with redshift) of ip. Using (2.29) one finds that the "equation of state parameter" 
w to t = -Ptot/ptot corresponding to the total energy and pressure (including now the kinetic 
contributions of (p; i.e. p tot = p + (m 2 P /2)(dip/dt) 2 , P tot = P + (m 2 P /2) (dp/di) 2 ) is given in 
terms of the "matter" equation-of-state parameter w = P/pbj 

X — uu 

w tot = w + ^-(cp') 2 ■ (2.30) 

The knowledge of u> to t then allows one to write explicitly the energy-balance equation dp to t + 
3(ptot + -Ptot)^l n o = 0, which is easily solved in the simple cases where w tot is (approximately) 
constant. 

We see from Eq. (2.28) that, during the radiation era (starting, say, immediately after the 
end of inflation), i.e. when the universe is dominated by an ultra- relativistic gas (pa~ 3Pa = 
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0), the "driving force" on the right-hand side of (2.28) vanishes, so that tp is not driven 
further away towards infinity. Actually, one should take into account both the "inertial" 
effect of the "velocity" tp' acquired during the preceding inflationary driving of tp, and the 
integrated effect of the many "mass thresholds", Ta ~ itla, when some component becomes 
non-relativistic (so that pa — 3Pa 7^ 0). Using the results of [4,3] one sees that, in our 
case, both these effects have only a small impact on the value of tp. Therefore, to a good 
approximation tp ~ tp CIlA until the end of radiation era. 

On the other hand, when the universe gets dominated by non-relativistic matter, one 
gets a non-zero driving force in Eq. (2.28). In the slow roll approximation, as the transient 
behaviour has died out, since w = P/p gets negligible, we have simply 

tp' m = -a m (tp) , (2.31) 

where tp' m stands the (^-velocity during matter domination, and a m (tp) denotes the coupling 
(2.27) to dark matter. 

The coupling to dark matter, a m (tp), depends on the assumption one makes about the 
asymptotic behaviour, at strong bare string coupling, of the mass of the WIMPs constituting 
the dark matter. One natural looking, minimal assumption is that dark matter, like all 
visible types of matter, is coupled in a way which levels off at strong-bare-coupling, as in 
Eq. (1.3). In other words, one generally expects that m m (tp) ~ m m (+oo)(l + 5 m e~ cip ) so that 
a m (tp) ~ — b m ce~ Cip . It is then easy to solve Eq. (2.31), with initial conditions tpo = yw, 
tp' = (inherited from radiation era) at the beginning of the matter era. But we shall not 
bother to write the explicit solution because it is easily seen that the smallness of e~ c,fcnd 
guarantees that the "driving force" oc a m (tp) remains always so small that the O(10) e-folds 
of matter era until vacuum-energy domination (or until the present) have only a fractionally 
negligible effect on tp. 

A more significant evolution of tp during the matter era is provided if, as first proposed in 
[11] and taken up in [25,26], dark matter couples much more strongly to tp than "ordinary" 
matter. Such a stronger coupling to dark matter, which is not constrained by usual equiva- 
lence principle experiments, follows assuming more general quantum corrections in the dark 
matter sector of the theory, i.e. corrections such that the dark matter mass m m (tp), instead 
of levelling off, either vanishes or keeps increasing at strong bare coupling: m m (<p) oc e CmV , 
so that a m = c m is a (negative or positive) constant. In [6] (but see also [12]) it has been 
shown that under the latter assumption (i.e. with a positive coupling parameter a m > 0) 
the dilaton can play the role of quintessence, leading to a late-time cosmology of acceler- 
ated expansion. By Eq. (2.31) we have tp = <^ en d — oi m p, where p is now counted from the 
end of the radiation era. Given that about nine e-folds separate us from the end of the 



13 



radiation era, we see that such an evolution might (if \a m \ is really of order unity) have a 
significant effect on the present value of tp (when compared with the value at the end of 
inflation, i.e. op cn d ~ ln(l /Sh) ~ 10). However, the running of tp during the matter era 
changes the standard recent cosmological picture and is therefore constrained by observa- 
tions. In fact, by equation (2.30), the total matter-era equation of state parameter iy tot in 
the presence of the dilaton reads w tot = (tp' m ) 2 /3. Accordingly, the matter density varies as 
p oc a~ 3( ~ 1+Wtot ^ = a-( 3 +(^m) 2 ) j possibly affecting the standard scenario of structure formation 
as well as the global temporal picture between now and the epoch of matter-radiation equal- 
ity. The compatibility with phenomenology therefore puts constraints on the magnitude of 
<p' m 2 = a m (tp) 2 . In [6] tutot < 0.1, i.e. v = <p' m 2 /0.3 < 1, was suggested to be the maximal 
deviation one can roughly tolerate during the matter era, the establishment of a more pre- 
cise bound being presently under study [13]. We shall therefore assume either that we are in 
the "normal" case where the dilaton does not couple more strongly to dark matter than to 
ordinary matter (so that a m ((p) ~ — b rn ce~ cip <C 1), or that, if it does, a 2 n < 0.3. This leads 
to a displacement of ip during matter era smaller than the dispersion tp end — y? c i,cnd ~ <£> c i,end 
produced by quantum fluctuations during inflation, Eq. (2.25). 

In this context one should also consider the attraction effect of a negative pressure com- 
ponent, either in the form of a (/^-dependent vacuum energy (dilatonic quintessence) or in 
the form of any other, (^-independent component (such as a "genuine" cosmological con- 
stant). Of course, the present recent (z ^ 1) accelerated expansion phase is very short (in 
"p-time") and sensible changes of the dilaton value since the end of matter domination are 
not expected. Still, it is crucial to estimate the present dilaton velocity tp' Q since it is related 
to the cosmological variations of the coupling constants (see next section). In the general 
case where both non-relativistic matter and (possibly (/^-dependent) vacuum energy density 
V((p) are present, the value of tp' predicted by our model is obtained by applying Eq. (2.28) 
(in the slow- roll approximation): 

(Q m + fV)(l - w )<p' = (^m + 2fV)y?o = -^m«m - 4:Q v a v . (2.32) 

In the above expression Q m and Qy are, respectively, the non-relativistic (dark) matter- 
and the vacuum- fraction of critical energy density (p c = {3/2)m 2 p H 2 ), and the already 
mentioned prescriptions ay = \d\nV(tp)/dtp, Py = —py = —V(ip) have been used. 

The value of ip' is therefore some combination of the values of a m and ay. We can have 
two classes of contrasting situations: In the first class, the dilaton couples "normally" (i.e. 
weakly) both to dark matter and to dark energy, i.e. both a m ~ — b m ce~ CLp 1 and ay 1 
and Eq. (2.32) implies <p' Q 1. In the second class, the dilaton couples more strongly to 
some type of dark matter or energy, i.e. either (or both) a m or/and ay is of order unity so 
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that ip' = 0(1). The second case is realized in the scenario of [6]. In the context of this 
scenario we have an exponential dependence of the potential on tp, V(<p) ~ V\e~ cip so that 
ay ~ — (c/4) and 

^ = cQ v - a m Q m < c ^ 
2f2y -|- f2 m 3 

The last inequality follows from the bound a m > c/2 (which is a necessary condition to have 
positive acceleration in the model [6]) and the reasonable bound fl m > 0.25. 

In the present work, we wish, however, to be as independent as possible from specific 
assumptions (as the ones used in [6]). Therefore, rather than insisting on specific (model- 
dependent) predictions for the present value of ip' we wish to find the (model-independent) 
upper bounds on the possible values of <p' set by current observational data. There are several 
ways of getting such phenomenological bounds, because the existence of a kinetic energy (and 
pressure) associated to dip/dt = Hip' has several observable consequences. A rather secure 
bound can be obtained by relating the value of ip' to the deceleration parameter q = —da/ a 2 . 
In the general class of models that we consider, the cosmological energy density and pressure 
have (currently) three significant contributions: dark matter ( Q m = p m /p c ), dark energy 
(fV) and the kinetic effect of a scalar field (f^ = pk/p c with p k = (m 2 P /2)(dip/dt) 2 = 
(m 2 P /2)H 2 p>' 2 so that Vt k = p' 2 /3). We assume (consistently with recent cosmic background 
data) that the space curvature is zero. Therefore we have the first relation 

tt m + tt v + tt k = 1 = tt m + tt v + <p' 2 /3 . (2.34) 

The deceleration parameter is given by the general expression 2q = T,a^a{^- + 3wa)- Using 
w m = 0, wy = —1 and w k = +1, we get 

4 

2q = Q m - 2Q V + -p' 2 . (2.35) 

Using the relation (2.34) above to eliminate £l v we get the following expression for ip' 2 in 
terms of the observable quantities q and Vt m 

p' 2 = l + q- hl m . (2.36) 

The supernovae la data [14] give a strict upper bound on the present value qo: qo < 0. 
A generous lower bound on the present value of Q m is Q m0 > 0.2 [15]. Inserting these two 
constraints in Eq.(2.36) finally yields the safe upper bound 

V?q 2 < 0.7, i.e. \ip' \ < 0.84 . (2.37) 
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To summarize, quite different rates of evolution for the dilaton are possible. A very 
slow variation is expected whenever dilaton couplings to both dark energy and dark matter 
follow the "normal" behaviour (1.3). Otherwise, dilaton variations on the Hubble scale 
are expected. However, cosmological observations set the strict upper bound (2.37) on the 
present time variation of tp. For the purpose of the present section (evaluating the current 
location of the dilaton) these two alternatives do not make much difference because the 
vacuum- dominance era has started less than about 0.7 e-folds away (ln(l + -2*) with < 1). 
Therefore, ip did not have enough "p-time" , during vacuum dominance, to move much, even 
if it is coupled to vacuum energy with ay ~ — (c/4) ~ 1. 

Finally, we conclude from this analysis that, to a good approximation (and using the 
fact that the phenomenology of the matter-era constrains the dark-matter couplings of the 
dilaton to be rather small), the value of ip now is essentially given by the value </? en d at the 
end of inflation, i.e. by Eq. (2.21). 

III. DEVIATIONS FROM GENERAL RELATIVITY INDUCED BY A RUNAWAY 

DILATON 

A. Composition-independent deviations from general relativity 

The previous section has reached the conclusion that present deviations from general 
relativity are given, to a good approximation, by the values of the matter-coupling coefficients 
q:a(p) given by Eq. (2.27) calculated at p ~ p cnd as given by Eq. (2.25). Let us now see the 
meaning of this result in terms of observable quantities. 

Let us first consider the (approximately) composition-independent deviations from gen- 
eral relativity, i.e. those that do not essentially depend on violations of the equivalence 
principle. Most composition-independent gravitational experiments (in the solar system or 
in binary pulsars) consider the long-range interaction between objects whose masses are es- 
sentially baryonic (the Sun, planets, neutron stars). As argued in [2,3] the relevant coupling 
coefficient a a is then approximately universal and given by the logarithmic derivative of the 
QCD confinement scale Aq CD ((/?), because the mass of hadrons is essentially given by a pure 
number times Aqcd (<£>)• [We shall consider below the small, non-universal, corrections to 
mA^p) and ola{p) linked to QED effects and quark masses.] Remembering from Eq. (1.1) 
the fact that, in the string frame (where there is a fixed cut-off linked to the string mass 
M s ~ {a')~ 1 / 2 ) the gauge coupling is dilaton-dependent (gp 2 = Bp ((f)), we see that (after 
conformal transformation) the Einstein-frame confinement scale has a dilaton-dependence of 
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the form 



AqcdM ~ C- 1 ' 2 B- l l\ V ) exp[-87r 2 63 1 B F (<p)] M s , (3.1) 

where 63 denotes the one-loop (rational) coefficient entering the Renormalization Group 
running of gv. Here B F ((p) denotes the coupling to the SU(3) gauge fields. For simplicity, 
we shall assume that (modulo rational coefficients) all gauge fields couple (near the string cut 
off) to the same B F (ip). This yields the following approximately universal dilaton coupling 
to hadronic matter 

ahadM^ In —2- +- -Z-W. (3.2) 



We recall that the quantity ahad(v 9 )> which measures the coupling of the dilaton to hadronic 
matter, should not be confused with any "strong" gauge coupling, a s = g 2 /^- Numerically, 
the coefficient in front of the R.H.S. of (3.2) is of order 40. Consistently with our basic 
assumption (1.3), we parametrize the f dependence of the gauge coupling g 2 F = Bp 1 as 

B~ F \ V ) =B F \+oo)[l-b F e-^]. (3.3) 



Note that, like b\ (see section 11 A), also b F is expected to be smallish [~ Bp 1 (+00) or, 
equivalently, ~ Cp 1 in the notations of (1.3)] and typically the ratio bp/b\ is of order unity. 
We finally obtain 

Ohadfa) -40 6 F ce" cv . (3.4) 



We can now insert the estimate (2.25) of the value of <p reached because of the cosmological 
evolution. Neglecting the 0{1) renormalization factor due to quantum noise, we get the 
estimate 

«had(^nd) -3.2-^5^, (3.5) 
0\ c 

aL d (^c„d)^10^ 2 ^. (3.6) 
As said above, it is plausible to expect that the quantity c (which is a ratio) and the ratio 

8 

bpjb\ are both of order unity. This then leads to the numerical estimate a^ ad ~ 10<5^- +2 , 
with 5h — 5 x 10~ 5 . An interesting aspect of this result is that the expected present value 
of a^ ad depends rather strongly on the value of the exponent n (which entered the inflaton 
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potential V(x) oc x n )- m the case n = 2 (i.e. V(x) = |m^x 2 ) we have a^ad ~ 2.5 x 10 8 , 
while if n — 4 (V(%) = | A % 4 ) we have aj; ad ~ 1.8 x 10~ 5 . 

How do these numbers compare to present (composition-independent) experimental limits 
on deviations from Einstein's theory [16]? This question has been addressed in the literature. 
Concerning solar-system (post-Newtonian) tests it was shown (see, e.g., [17]) that the two 
main "Eddington" parameters 7 — 1 and (3 — 1 measuring post-Newtonian deviations from 
general relativity are linked as follows to the dilaton coupling ahad(v 9 ) : 

l-l = -2-^-^-2oi td , (3.7) 
1 + (*La 



a -1 1 a had a had ^ 1 „ / „ 2 /, o\ 

^ ~ 2 (iT^Ld^ 2 had had ' (3 ' 8) 

where a' had = da had ((p)/d(p. 

From Eq. (3.4) we see that a' had ~ — c« had , so that the deviation f3 — 1 is 0(0:^) 
and thereby predicted to be too small to be phenomenologically interesting. This leaves 
7 — 1 ~ —2 «j 2 iad as the leading observable deviation. The best current solar-system limit on 
7 — 1 comes from Very Long Baseline Interferometry measurements of the deflection of radio 
waves by the Sun and is (approximately) (7 — 1| % 2x 1CT 4 , corresponding to a^ ad < 1CT 4 (see 
[16] for reviews and references). In addition to solar-system tests, we should also consider 
binary-pulsar tests which provide another high-precision window on possible deviations from 
general relativity. They have been analyzed in terms of the two quantities ahad (denoted a) 
and a had (denoted 0) in [18]. The final conclusion is that the binary-pulsar limit on a had is 
of order a^ ad < 1CT 3 . 

At this stage it seems that the runaway scenario explored here is leading to deviations 
from general relativity which are much smaller than present experimental limits. However, 
we must turn our attention to composition- dependent effects which turn out to be much more 
sensitive tests. 



B. Composition-dependent deviations from general relativity 

Let us then consider situations where the non-universal couplings of the dilaton in- 
duce (apparent) violations of the equivalence principle. Let us start by considering the 
composition-dependence of the dilaton coupling a a, Eq. (2.27), i.e. the dependence of a a on 
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the type of matter we consider. The definition of a a is such that, at the Newtonian approx- 
imation, the interaction potential between particle A and particle B is — Gab ?tia ttib / r ab 
where [3] 

Gab = G{l + a A a B ). (3.9) 



Here, G is the bare gravitational coupling constant entering the Einstein-frame action (2.4), 
and cxa = «a(v 9 ) is the strength of the dilaton coupling to A-particles, taken at the present 
(cosmo logically determined) VEV of (p. The term a a c*b comes from the additional attractive 
effect of dilaton exchange. Two test masses, made respectively of A- and B-type particles 
will then fall in the gravitational field generated by an external mass m E with accelerations 
differing by 

fAa\ n a A -a B 

— =2 ■ ~ {a A - ol b ) a E ■ (3.10) 

\ a / ab a,A + a B 

We have seen above that in lowest approximation a a — «had does not depend on the com- 
position of A. We need, however, now to retain the small composition-dependent effects 
to a a linked to the (^-dependence of QED and quark contributions to itia- This has been 
investigated in [3] with the result 



(tL-(^)>M#) +c WsW( 



-) 

mJ 



AB 



(3.11) 



where (AX)ab = Xa — X B , where B = N + Z is the baryon number, D = N — Z the 
neutron excess, E = Z(Z — 1)/ (N + Z) 1 ^ 3 a quantity linked to nuclear Coulomb effects, and 
where M = m/u denotes the mass in atomic mass unit, u = 931.49432 MeV. It is difficult 
(and model-dependent) to try to estimate the coefficients C B and C B . It was argued in [3] 
that their contributions to (3.11) is generically expected to be sub-dominant with respect 
to the last contribution, oc C B , which can be better estimated because it is linked to the 
(^-dependence of the fine-structure constant e 2 oc Bp 1 (ip). This then leads to the numerical 
estimate C E — 3.14 x 1CT 2 and a violation of the universality of free fall approximately given 
by 



V a J ab 



E\ f E 
m)a~\M 



(3.12) 



The values of B/M, D/M and E/M have been computed in [19]. For mass-pairs that 
have been actually used in recent experiments (such as Beryllium and Copper), as well as 
for mass-pairs that are planned to be used in forthcoming experiments (such as Platinum 
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and Titanium) one finds: (E/M) Cu - {E/M) Bc = 2.56, {E/M) Pt - {E/M) Ti = 2.65. Using 
the average estimate A(E/M) ~ 2.6, we get from (3.12) and (3.6) the estimate 



(^) ~ 5.2 x 10- 5 < d ~ 5.2 x 10- 4 (J^J 5^ . (3.13) 

Note also (from (3.7)) the link between composition-dependent effects and post-Newtonian 
ones 

— ) ~ -2.6 x lCr 5 ( 7 - 1). (3.14) 

As current tests of the universality of free fall (UFF) have put limits in the 10~ 12 range (e.g. 
(Aa/a) Be cu = (—1.9 ± 2.5) x 1(T 12 from [20]), we see from Eq. (3.14) that this corresponds 
to limits on 7 — 1 or a 2 ad in the 10~ 7 range. Therefore tests of the UFF put much more 
stringent limits on dilaton models than solar-system or binary-pulsar tests. 

If we insert the estimate 5h ~ 5 x 10~ 5 in (3.13) we obtain a level of violation of UFF 
due to a runaway dilaton which is 

A " 1.3 f-^Vx 10- 12 forn = 2, (3.15) 



a \b\c 



— ~ 0.98 ( p-\ x 10" 9 for n = 4 . (3.16) 
a \b x cj 

At face value, one is tempted to conclude that a scenario with n = 4 (i.e. V(x) oc x 4 ) tends 
to be too weak an attractor towards tp = +oo to be naturally compatible with equivalence- 
principle tests. [See, however, the discussion below.] On the other hand, the simple scenario 
n = 2 (V(x) = \ ™ 2 X X 2 ) is quite appealing in that it naturally provides enough attraction 
towards tp = +oo to be compatible with all existing experimental tests. At the same time it 
suggests that a modest improvement in the precision of UFF experiments might discover a 
violation caused by a runaway dilaton. 



C. Cosmological variation of "constants" 

Let us now consider another possible deviation from General Relativity and the standard 
model: a possible variation of the coupling constants, most notably of the fine structure 
constant e 2 /he on which the strongest limits are available. We will discuss first the effects 
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due to the cosmological time-variation of the homogeneous component of ip and, in the next 
subsection, the possible spatial (and time) variations due to quantum fluctuations of tp as 
they got amplified during inflation. 

Consistently with our previous assumptions we expect e 2 oc Bp 1 ^) so that, from (3.3), 

e 2 (^) = e 2 (+oo)[l-& F e-^]. (3.17) 

The present logarithmic variation of e 2 (using again dp = H dt; ip' = dip /dp) is thus given by 

d In e 2 d In e 2 



Hdt ~~ dp 



~b F ce~ Cip <p' , (3.18) 



where the current value of ip', ip' , is given in general by Eq. (2.32). Using Eq. (3.4), we can 
rewrite the result (3.18) in terms of the hadronic coupling: 

dine 2 1 , /.-, ,«\ 

W -40^o- (3-19) 

As said in section II B, we have basically two alternatives concerning the current coupling 
of the dilaton to the dominant energy sources in the universe. These two alternatives lead 
to drastically different predictions for the current value of the rate of variation of the fine- 
structure constant. We shall consider these two alternatives in turn. 

In the conservative case where the dilaton does not play any special role in the present 
accelerated phase of the universe (ay — 0) nor does it have any stronger coupling to dark 
matter than to visible matter (a m ~ —b m ce~ C(p ) the dilaton "velocity" ip' is exponentially 
suppressed ( so that, from (2.34), fl v ~ 1 — Q m ) and by Eq. (2.32) one obtains 

' r ' [l! ~ b F ce-^a m {ip) ~ -^-b F b m c 2 e~ 2c *. (3.20) 



Hdt fl m + 2Q v r mvr/ 2-Q 



in 



An indicative value for the ratio VL m /(VL m + 2VL V ) ~ fi m /(2 — f2 m ), by taking for instance 
VL m = 0.3, is 0.18. As above, it is useful to relate (3.20) to the estimate (3.4) for «had- This 
yields 

dlne2 1 Q rn b m , 

Hdt ~ (40) 2 2 -n m b F a ^ d - [6 ] 

In terms of the UFF level Aa/a predicted by our model in (3.13) we see also that 

Hdt 2 - n m b F a { ' ' 
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Even if the universe were completely dominated by dark matter (fi m = 1) we see, assuming 
that b m /bp is of order unity, that current experimental limits on UFF (Act/a < 10 -12 ) imply 
(within dilaton models) that \d\ne 2 /dt\ < 10~ n H ~ 10~ 21 yr _1 (the sign of d In e 2 / dt being 
given by the sign of b m jbp). This level of variation is much smaller than the current best 
limit on the time variation of e 2 , namely \dlne 2 /dt\ < 5 x 10 _17 yr _1 ~ 5 x 10~ 7 H, as 
obtained from an analysis of Oklo data [21]. (Note that the assumption-dependent analysis 
of Ref. [22] gives a limit on the variation of e 2 which is strengthened by about two orders of 
magnitude.) 

The situation, however, is drastically different if we consider the alternative case where 
the dilaton coupling to the current dominant energy sources does not tend to triviality, as 
in the case of a 92-dependent vacuum energy V(tp) — Vq + V\e~ Cif when the first term is zero 
or negligible. In such a case the dilaton shares a relevant part of the total energy density 
and more significant (though still quite constrained by UFF data) variations of the coupling 
constants are generally expected. A general expression for the dilaton "velocity" is given in 
eq. (2.36) in terms of observable quantities. Using Eqs. (2.36) and (3.19) one can relate the 
expected variation of the electromagnetic coupling constant to the hadronic coupling: 



We can also use the estimate (3.5) relating a ha d to the density fluctuations generated 
during inflation. We obtain 



However, in view of the theoretical uncertainties attached to the initial conditions Xm 
and (pi n used in the estimate (3.5), as well as the ones associated to the order unity ratio 
bp/ipxc), it is more interesting to rewrite our prediction in terms of observable quantities. 
Using again the link Eq. (3.13) between a had and the observable violation of the universality 
of free fall (UFF) the above result can be written in the form 



Note that the sign of the variation of e 2 is in general model-dependent (as it depends 
both on the sign of bp and the sign of (p' ). Specific classes of models might, however, favour 
particular signs of de 2 /dt. For instance, from the point of view of [5] one would expect the 
0(e~^) terms in Eq. (1.3) to be positive, which would then imply that bp is positive. If we 




(3.23) 




(3.24) 




(3.25) 
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combine this information with the prediction Eq. (2.33) of the model [6] implying that tp' is 
also positive, we would reach the conclusion that e 2 must be currently increasing. 

Independently of this question of the sign, we see that Eq. (3.25) predicts an interesting 
link between the observational violation of the UFF (constrained to Aa/a ^ 1CT 12 ), and 
the current time- variation of the fine-structure constant. Contrary to the relation (3.22), 
obtained above under the alternative assumption about the dilaton dependence of the domi- 
nant cosmological energy, which predicted a relation linear in Aa/a, we have here a relation 
involving the square root of the UFF violation (such a relation is similar to the result of [11] 
which concerned the time- variation of the Newton constant). 

The phenomenologically interesting consequence of Eq. (3.25) is to predict a time- 
variation of constants which may be large enough to be detected by high-precision labora- 
tory experiments. Indeed, using Hq ~ 66 km/s/Mpc, and the plausible estimates Q m = 0.3, 
q = -0.4, Eq. (3.25) yields the numerical estimate d In e 2 /dt ~ ±0.9xl0 -16 ^lO^Aa/ayr -1 . 
Therefore, the current bound on UFF violations (Aa/a ~ 10~ 12 ) corresponds to the level 
10 _16 yr _1 , which is comparable to the planned sensitivity of currently developed cold-atom 
clocks [28]. [Present laboratory bounds are at the 10~ 14 yr _1 level [29,28].] Note that if 
we insert in Eq. (3.25) the secure bounds Q m > 0.2 and go < (leading to the limit 
Eq. (2.37)), we get as maximal estimate of the time variation of the fine-structure constant 
d\iae 2 /dt ~ ±2.0 x 10~ 16 y / 10 12 Aa/ayr~ 1 . We note also that the upper limit on the varia- 
tion of e 2 given by the Oklo data, i.e. |<ilne 2 /<it| < 5 x 10~ 17 yr _1 [21], "corresponds" to a 
violation of the UFF at the level ~ 10~ 13 . 

In this respect, it is interesting to consider not only the present variation of e 2 (the 
only one relevant for laboratory experiments), but also its variation over several billions of 
years. (We recall that the Oklo phenomenon took place about two billion years ago, and 
that astronomical observations constrain the variation of e 2 over the last ten billion years 
or so). In particular, an interesting question is to see whether our model could reconcile 
the Oklo limit (which corresponds to a redshift z ~ 0.14) with the recent claim [24] of 
a variation Ae 2 /e 2 = (—0.72 ± 0.18) x 10 -5 around redshifts z ~ 0.5 — 3.5 as proposed 
in [25,26]. The only hope of reconciling the two results would be to allow for a faster 
variation of e 2 for redshifts z > 0.5. Such recent redshifts have (apparently) been connected 
to a transition from matter dominance to vacuum dominance. Let us see whether taking 
into account this transition might allow for a large enough change of e 2 around redshifts 
z w 0.5 — 3.5. We must clearly assume the "strong coupling" scenario a m — 0(1). In 
this scenario, the variation of y? during the matter era is given by Eq.(2.31). Neglecting, 
for simplicity, the transient evolution effects localized around the matter-vacuum transition 
(and treating both ip' m = —a m and ip' v = ip' Q as constants), the solution giving the recent 
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cosmological evolution of ip reads ip — </?o = — v^o l^Cl + z ) during the vacuum era, and 
(p — (p = — (/?q ln(l + z*) — ip' m ln[(l + z)/(l + z*)] during the matter era (the index refers 
to the present epoch, i.e. z = 0; z* denotes the transition redshift). Inserting this change in 
Eq.(3.17) leads to the following expression for the cosmological change of the fine-structure 
constant: 



^ = -sign(o F ) 3.5 x 1(T 6 [<p' ln(l + z.) + p>' m In -^—\ WlO 12 — . (3.26) 

6q 1 + Z* V 0, 

Here, we have written the result for the matter era. During the vacuum era the bracket 
is simply [(/9' ln(l + z)\. Remembering that the absolute value of (p' m is (like that of (p' v ) 
observationally constrained to be smaller than v^3 ~ 0.55 (and that (f' is also constrained 
by \<Pv\ < 0.84), we see that there is no way, within our model, to explain a variation of 
e 2 as large as Ae 2 /e 2 = (-0.72 ± 0.18) x 10~ 5 around redshifts z w 0.5 - 3.5 [24]. In 
our model, even under the assumption that UFF is violated just below the currently tested 
level, such a change would have to correspond to a value \ip' m \ > 2, entailing observationally 
unacceptable modifications of standard cosmology. [For instance, in the model [6] a value 
as large as a m > 1 already leads to a pathological behaviour ("total dragging") where all 
the components scale like radiation.] This difficulty of reconciling the Oklo limit with the 
claim of [24] was addressed in [26,25] within a different class of models, namely with a field 
4> which does not couple universally to all gauge fields F^ v , as the dilaton ip is expected 
to do. The fact that the field in [26] (or ip in [25]) is assumed to couple only to the 
electromagnetic gauge field drastically changes our Eq. (3.13) and allows one to satisfy the 
UFF limit Aa/a < 10~ 12 , for a stronger coupling of <fi to electromagnetism than in our 
class of models, i.e. (in our notation) for a larger dhi B p{(p) / dip . This explains why Ref. 
[26] could construct some explicit (but fine-tuned) models in which all observational limits 
(UFF, Oklo,...) could be met and still allow for a variation of e 2 as strong as the claim[23]. 
The maximal variation predicted by Eq.(3.26) for redshifts corresponding to the matter era 
(obtained when Aa/a = 10~ 12 and ip' m = ±V0.3; and assuming a smaller value of ip' to 
be compatible with the Oklo constraint), is of order Ae 2 /e 2 = ±1.9 x 10~ 6 . This is only a 
factor ~ 4 below the claim [24] and is at the level of their one sigma error bar. Therefore a 
modest improvement in the observational precision (accompanied by an improved control of 
systematics) will start to probe a domain of variation of constants which, according to our 
scenario, corresponds to an UFF violation smaller than the 10~ 12 level. 
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D. Spatio-temporal fluctuations of the "constants" 



We now turn to the second possible source of spatial-temporal variations for e 2 in our 
model, the quantum fluctuations of the dilaton generated during inflation. Within linear 
perturbation theory, the relevant calculation may be summarized as follows. 

Consider a flat FRW universe ds 2 = —dt 2 + a(t) 2 T,idx 2 . The dilaton fluctuations can be 
expanded in Fourier components 5ipk of given comoving momentum k as follows: 

6(p(x,t) = - * / d 3 k5 Vk (t)e lk *, (3.27) 



(2tt) 3 /2 

where t is the cosmological time. Each Fourier mode 5^ "leaves" the horizon during inflation 
with an amplitude ~ H cx (k) /y/2k 3 [27] where, by definition, H ex (k) is the value of the 
dimensionless Hubble expansion rate as ka~ l equals H during inflation (note that we denote 
here H ex what was denoted H x above). Well after the exit (k <C aH) the amplitude of each 
mode "freezes out", i.e. remains roughly constant, until it reenters the horizon during the 
post-inflationary epoch {ka~^ ~ H re ). After re-entry the amplitude starts to damp out as 
a -1 . For a given Fourier mode 5ipk(t), the latter damping effect is described by the piecewise 
function 

' 1 if a^H^k < (z + l) 1 ' 2 

f z (k) = J a 2 H 2 (z + l)k~ 2 if (z + If/ 2 < cl^Hq 1 k < 10 2 (3.28) 

k 10~ 2 a H (z + l)^- 1 if a^H^k > 10 2 

Here the cosmological redshift z = ao/a(t) — 1 has been introduced in replacement of the 
cosmological time t. The first case refers to Fourier modes that have not reentered yet at 
redshift z and whose amplitudes are still frozen. The second and third cases refer to modes 
that reenter during matter and radiation domination respectively. Putting all together, and 
assuming a gaussian probability distribution for the perturbations, we have: 

H 2 (k) 

(5^(ty 5<pv(t')) = -sii /,(*) f Z '(k) 5 3 (k - k') . (3.29) 

Possible spatial/temporal variations of e 2 induced by the fluctuations of the dilaton will be 
given by 



e 2 



= ^A fl >| (x , w) , (3.30) 

(x,t;x',t') ai P 



where the r. m.s A flu V 

between two events (x, t) and (x', t') is defined as follows 
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I 2^ [ /z(A;)2 + fAk)2 ~ 2f *^ f *'W 



(27r) 3 7 2 k 3 
2 I 



1 - 



sin /err 



/ex 



(3.31) 
(3.32) 



Here, x = |x — x'| is the coordinate distance between the two events and, consistently with 
the slow-roll approximation, the Hubble expansion rate at exit has been assumed to be 
scale- invariant: H cx {k) ~ H cx ~3x 1CT 5 . 

If one considers spatial fluctuations over terrestrial or solar system proper length scales 
I = aok' 1 <C Hq 1 at the present time t — t' — t , the first square brackets in (3.32) vanishes 
and one can expand the sine function at small kx obtaining 

H cx HqI 



A flue I 
f\l;z=0 — 



A fluc e 2 



2tt 



~ 10 a had H cx H l. 



(3.33) 



l;z=0 



As expected, these variations are extremely small, A fluc e 2 /e 2 |z ;z= o — 10 _33 //km. It is also 
interesting to compare dilaton fluctuations at different redshifts along a comoving observer 
wordline. By putting x ~ in (3.32) the second term in the square brackets vanishes and 
one has: 

1/2 



A flue I ^ -^CX 1 

A (p\ z;x=0 ~ 



2tt s/2 



log(l + z) - 



l + z 



+ 



10" 



2-J + 



(3.34) 



It is slightly more complicated to compare dilaton fluctuations between "now" and events 
at redshift z along a null ray. Expanding in powers of z around z = one gets from (3.32), 
after a straightforward calculation: 



flue. 



H cx 

2vr 



(3.35) 



1, the effects of dilatonic fluctuations give A (p\ z =i 



. flue, 



Numerically, at redshift z 
H ex /(2ir) ~5x 10 -6 . This is to be contrasted with the effects of the cosmic, homogeneous 



evolution which yields Atp\ 



2 = 1 



ot r 



In the "normal" case where a r 



-C(f 



S 



4/(rt+2) 
H > 



the two effects, though a priori unrelated, are related in our scenario , when n — 2. In- 
deed, if n = 2, <5^ /(n+2) = S H ~ 5 x 10~ 5 is linked to H cx /(2n) via 5(Xex) = AH cx /(2n) 
with A = (8/3)V/d x V = (8/3) (x/n) ~ 40/(301) ~ 10. On the other hand, in the case 
where tp is strongly coupled to dark matter, the homogeneous evolution Atp\ z= i ~ a m ~ 1 
is parametrically larger than the fluctuations A fluc y?| 2= i ~ H cx /(2n). 
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To conclude on this subsection, we see that the inhomogeneous space-time fluctuations 
of the fine-structure constant are typically too small to be observable (if the limits from 
UFF are already satisfied), being suppressed, relative to their natural values HqI , Hot, by 
the small factor ah&dH cx . 



IV. SUMMARY AND CONCLUSION 



We have studied the dilaton-fixing mechanism of [3] within the context where the dilaton- 
dependent low-energy couplings are extremized at ip = +00, i.e. for infinitely large values 
of the bare string coupling g 2 s = ~ e ctp . [The crucial coupling to the inflaton, say \(<p) 
in Eq.(2.11), must be minimized at ip — > +00; the other couplings can be either minimized 
or maximized there.] This possibility of a fixed point at infinity (in bare string coupling 
space) has been recently suggested [5], and its late-cosmological consequences have been 
explored in [6]. We found that a primordial inflationary stage, with inflaton potential V(x) = 
X(ip) X n /n, was much less efficient in decoupling a dilaton with least couplings at infinity 
than in the case where the least couplings are reached at a finite value of (p (as in [3,8]). 
This reduced efficiency has interesting phenomenological consequences. Indeed, it predicts 
much larger observable deviations from general relativity. In the case of the simplest chaotic 
potential [10] V(x) = \ m x( { P)x 2 ^ we find that, under the simplest assumptions about the 
pre-inflationary state, this scenario predicts violations of the universality of free fall (UFF) 
of order Aa/a ~ 5 x 10 -4 <5^ where 5h is the density fluctuation generated by inflation on 
horizon scales. The observed level of large-scale density (and cosmic microwave background 
temperature) fluctuations fixes Sh to be around 5 x 10~ 5 which finally leads to a prediction for 
a violation of the UFF near the Aa/a ~ 10~ 12 level. This is naturally compatible with present 
experimental tests of the equivalence principle, and suggests that a modest improvement in 
the precision of UFF tests might be able to detect a deviation linked to dilaton exchange 
with a coupling reduced by the attraction towards the fixed point at infinity. Because of the 
presence of unknown dimensionless ratios (c, bp/b\) in our estimates, and of quantum noise in 
the evolution of the dilaton, we cannot give sharp quantitative estimates of Aa/a. However, 
we note that dilaton-induced violations of the UFF have a rather precise signature with a 
composition-dependence of the form (3.11), with probable domination by the last (Coulomb 
energy) term [3]. As explored in [19] this signature is quite distinct from UFF violations 
induced by other fields, such as a vector field. We note that the approved Centre National 
d'Etudes Spatiales (CNES) mission MICROSCOPE [30] (to fly in 2004) will explore the 
level Aa/a ~ 10~ 15 , while the planned National Aeronautics and Space Agency (NASA) and 
European Space Agency (ESA) mission STEP (Satellite Test of the Equivalence Principle) 
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[31] could explore the Aa/a ~ 10~ 18 level. Our scenario gives additional motivation for 
such experiments and suggests that they might find a rather strong violation signal, whose 
composition-dependence might then be studied in detail to compare it with Eq. (3.11). 

In the case of inflationary potentials V(x) oc x™ with n > 2 our simplest estimates predict 

8 

a violation of the UFF of order Aa/a ~5x lCr 4 5^ +2 which is larger than 1CT 12 . At face 
value this suggests that existing UFF experimental data can be interpreted as favouring 
n < 2 over n > 2. However, we must remember that our estimates have made several 
simplifying assumptions. It is possible that the large quantum fluctuations of the inflaton in 
the self- regenerating regime x > Xm, with Xm defined by Eq. (2.18), can give more time for tp 
to run away towards large values, so that the effective value of e C(pin to be used in Eq. (2.21) 
turn out to dominate the first term in the R.H.S. that we have used for our estimates. We 
leave to future work a study of the system of Langevin equations describing the coupled 
fluctuations of and x during the self-regenerating regime. 

Finally let us note some other conclusions of our work. 

We recover the conclusion of previous works on dilaton models that the most interest- 
ing experimental probes of a massless weakly coupled dilaton are tests of the UFF. The 
composition-independent gravitational tests (solar-system, binary-pulsar) tend to be much 
less sensitive probes (as highlighted by the relations (3.14), (3.21) and (3.22)). 

However, a possible exception concerns the time-variation of the coupling constants. 
Here the conclusion depends crucially on the assumptions made about the couplings of the 
dilaton to the cosmologically dominant forms of energy (dark matter and/or dark energy). 
If these couplings are of order unity (and as large as is phenomenologically acceptable, i.e. 
so that ((f'o) 2 = 0.7), the present time variation of the fine-structure constant is linked to the 
violation of the UFF by the relation d\ne 2 /dt ~ ±2.0 x 10~ 16 ^lO^Aa/ayr" 1 . [The most 
natural sign here being +, i.e. bp > 0, which corresponds to smaller e 2 in the past, just 
as suggested by the claim [24].] Such a time variation might be observable ( if Aa/a is not 
very much below its present upper bound ~ 10~ 12 ) through the comparison of high-accuracy 
cold-atom clocks and/or via improved measurements of astronomical spectra. 

More theoretical work is needed to justify the basic assumption (1.3) of our scenario. 
In particular, it is crucial to investigate whether it is natural to expect that the sign of the 
crucial coefficient b\ in Eq. (2.12) be indeed positive. [Recall that the general mechanism of [3] 
is an attraction towards "Least Couplings" while Eq. (1.3) with 0(e^ (/> ) > leads to largest 
couplings at infinity] Note in this respect that the sign of the other b^s is not important 
as, once inflation has pushed e CLp to very large values e ciPcnd , the subsequent cosmological 
evolutions tend to be ineffective in further displacing <p. 
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APPENDIX A: THE STOCHASTIC EVOLUTION OF THE DILATON 



In this appendix we study the stochastic evolution of the dilaton ip during inflation as 
described by the Langevin-type equation (2.23). We restrict our attention to the region of 
phase space where the evolution of the inflaton x is classical, and to a power-law potential 
of the form (2.11). It follows that the inflaton evolves according to the classical slow-roll 
equation (2.14) whose solution reads 

x 2 = xL - n P> ( A1 ) 



where p, the parameter defined in (2.7), is shifted in such a way that p in = 0. Equation 
(2.23) takes the form 

| = ^Ace^+e(p), (A2) 

where is a gaussian stochastic variable (GSV), with a "time-dependent" r.m.s. amplitude 
H(p)/2n: 

<£(PiK(P2)> = (|}2*(Pi-P3), ( A3 ) 

[the relation with the normalized random white noise term of (2.23) is = ^2(p)H /2%]. 
For any given source term the formal solution of (A2) reads 



e c V ( P ) = e c Vin e c V ( P ) + h _><? r p dp > e c[ V (p)-v( P ')] ; = r 'dp' t(p') . (A4) 

2 Jo jo 



Note that the classical solution in (2.17), e CLPcl( - p ^ = e CLPin + (b\c 2 /2) p , can be easily recovered 
in the small noise limit £(p) — > 0, i](p) — > 0. 
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It proves convenient to compare the true solution to the classical one by studying the 
statistical behaviour of the ratio A(p) = e CLp<yP ^ /e cificl ^ p \ As we will show below, (e cv ^) = 
(9(1). Moreover, we are also assuming e c<fin = 0(1) or, at least, e ctp ™ <C (b\c 2 /2)p (see 
Section II for details) so that the leading contribution to the first equation in (A4) is given 
by the integral, and we have 



A(p) 



- [ dp'e 
p Jo 



>f,c[v(p)-v(p')] 



(A5) 



Since £(p) is a GSV, also its integral r](p) is a (centered) GSV. Moreover, if x is a GSV 
with a 2 = (x 2 ) — (x) 2 = (x 2 ), by Bloch's theorem y = e x is a new stochastic variable with 
(y) = (e x ) = e^l 2 and a 2 = e 2 ^ — . The average value of A(p) thus reads 



1 fP 

(A(p)) ~ - / dp e 
p Jo 



l (c 2 /2) (WP)~V(P')} 2 ) 



(A6) 



The exponent on the right hand side of the above equation can be estimated by using (A3) 
and the slow-roll approximation H 2 ~ 2V(x, V 9 )/3 = 2A(</?)x n /3n ~ 2\ OQ x n /3n. One gets 



([flip) - vtf)Y) = 



(2vr) 2 



H 2 dp" 



n 



2(n + 2) 



XVP ) \ / XVP) \ 



Xi 



(AT) 



where Xin is the value at exit from self-regenerating inflation: H(xm)/^ = n/(2xm) (see 
Section II for more details). Since we are interested in evaluating (A6) at the end of inflation, 
P = Pcnd — x\J n i we can thus write 



{[riip) - v(p')V) 



n 



2(n + 2) 



1 



p 

Pend 



l (n+2)/2 



(A8) 



When evaluated at p' = 0, the above formula gives (r](p) 2 } = n/(2(n + 2)). Thus the 
normalization factor to the initial condition in (A4) is of order one, as anticipated: (e CT/ ^) = 
e i c Mp) 2 } = o(\y prom (A6) and (A8) we have: 



(A(p end )) 



1 fPend 



Pend 

1 



/•Pond 

/ dp' 
Jo 



cxp 



c n 



p 



(n+2)/2' 



4(n + 2) V p end/ 







exp 



c z n 



4(n + 2) 
c 2 n9 \ 



x 



(n+2)/2 



dx 



(A9) 



exp 



^4(n + 2) 



with < 9 < 1. 
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We can estimate the dispersion of the same quantity by expanding the exponential inside 
the integral (A6) in powers of 

A(p) ~ 1 + J^d^(jf) + ~ J" dp' ^ dp'W')) 2 + ■ ■ ■ ■ (A10) 

At lowest order in the variance of the above quantity calculated at p 

^>=(gf*'|V^)) 2 ) 
= 4/V/V/* 

p z Jo JO Jmax(p',p") (27TJ 2 

p 2 Jo Jp> (27TJ 2 

As in equation (A7) we can use the slow-roll approximation and obtain 

= ^ 1 1 f 1 - *> , " +2,/2 < fa = c2 °< 1 ) ■ < A12 > 



= Pend reads 



(All) 
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